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I .  Introduction 

The  Naval  Postgraduate  School  has  started  a  project  to 
design  and  build  a  mini-satellite  (ORION)  .  This  mini-satellite 
bus  will  support  a  small  payload  for  various  space  science, 
technology  and  military  missions.  The  first  specific  payload  is 
not  determined,  yet,  and  therefore,  its  major  design  parameters, 
such  as,  launch  vehicle,  attitude  stabilization  scheme,  are  not 
known.  A  feasibility  study  and  a  preliminary  design  of  the  mini- 
satellite  were  performed  by  Boyd  (Ref .  1)  .  In  this  preliminary 
design  the  launch  of  the  satellite  from  the  GAS  (Get-Away- 
Specials)  canister  of  a  Space  Shuttle  and  single  spin 
stabilization  were  assumed 

This  preliminary  design,  as  defined  in  Ref.  1,  is 
investigated  from  the  structural  dynamics  and  attitude  dynamics 
points  of  view.  These  include  the  survivability  of  the  satellite 
under  severe  shock  and  vibration  environment  during  launch 
periods,  and  the  stability  of  the  single  spin  attitude 
stabilization  of  the  satellite.  A  rigorous  investigation  of  the 
survivability  satellite  under  shock  and  vibration  environment  is 
not  possible  in  this  preliminary  design  stage  since  the 
acceptible  acceleration  level  of  each  electronic  and  mechanical 
components  and  detail  design  of  connections  and  fasteners  are  not 
defined.  However,  it  seems  that  the  present  design  of  bolts  and 
nuts  connection  and  no  damping  treatment  (Figure  1)  contains  high 
failure  probability  areas  considering  the  wide  input  frequency 
range  (20-2000  Hz)  of  random  vibration  in  GAS  Canisters  (Figure 
2).  " 

The  main  body  of  the  NPS  mini-satellite  is  long  cylindrical 
shape  (See  Figure  1) .  The  principal  moment  of  inertia  about  the 
axial  axis  of  symmetry  is  smaller  than  the  principal  moments  of 
inertia  about  the  radial  axes  of  the  main  body  without  the  four 
radial  booms.  It  is  well  known  that  the  attitude  of  a  rigid  body 
spinning  about  the  axis  of  the  minimum  principal  moment  of 
inertia  is  unstable.  Therefore,  in  the  preliminary  design  of  the 
NPS  satellite,  the  lengths  of  booms  and  the  masse  of  the 
magnetometer  at  the  end  of  each  boom  are  such  that  the  total 
moment  of  inertia  about  the  axial  axis  is  the  maximum,  assuming 
the  four  radial  booms  are  rigid.  However,  long  booms  are 
flexible  and  the  flexibility  of  booms  affects  the  attitude 
stability  of  a  spinning  body  negatively.  Therefore,  each 
specific  design  of  a  spin  stabilized  satellite  with  flexible 
parts  should  be  investigated  to  determine  the  stability 
conditions   between  design  parameters. 

Dynamic  equations  governing  3-dimensional  motions  of  a  body 
(especially  a  flexible  body)  in  space  are  highly  nonlinear  and 
complicated.       Therefore,     a    free    or    forced    motion    and    the 


stability  of  the  motion  of  a  body  in  3-dimensional  space  can  be 
best  investigated  by  the  computer  simulation  for  each  specific 
case  using  a  highly  capable  computer  code  such  as  DISCOS  (Ref. 
2)  ,  or  TREETOPS  (Ref.  3) .  In  this  preliminary  design  stage,  and 
without  an  available  code,  the  stability  of  the  simple  spinning 
motion  of  the  NPS  mini-satellite  about  the  axial  axis  of  symmetry 
is  investigated  analytically  in  Chapter  II.  The  analytical  study 
of  the  spin  stability  of  free  motions  of  elastic  disspative 
systems  was  performed  previously  by  Teixeira-Filho  and  Kane  (Ref. 
4).  The  results  of  this  study  were  applied  to  the  stability 
studies  of  a  spinning  rigid  satellite  that  carries  four 
elastically  mounted  rigid  antennas  (Ref.  5)  and  of  a  spacecraft 
consisting  of  a  rigid  body  and  four  elastic  booms  using  finite 
number  of  beam  elements  (Ref.  6)  .  In  the  present  investigation 
of  the  stability  of  free  spinning  motion  of  the  NPS  mini- 
satellite  the  basic  results  of  the  Ref.  4  are  applied,  using 
natural  vibration  modes  of  the  booms  as  the  admissble  functions 
of  internal  coordinates. 


II.    Analysis 

II-l   system  Description 

The  NPS  mini-satellite  S  to  be  analyzed  consists  of  a  rigid 
body  B  of  cylindrical  shape  and  four  flexible  booms  C,  D,  E,  F, 
as  shown  in  Figure  3.   B*  is  the  mass  center  of  B  and  X-  ,  X_ ,  X_ 

are  the  principal  axes  of  inertia  of  B.   When  undeformed,  boom  C 
and  D  lie  on  the  opposite  side  of  B  in  the  X  -X-  plane  and  are 

parallel  to  X  while  booms  E  and  F  lie  on  opposite  sides  of  B  in 

the  X  -X_  plane  and  are  parallel  to  X_  .   The  four  booms  are 

assumed  to  be  rigidly  attached  on  the  cyclindrical  surface  of 
radius  R  of  body  B  in  the  radial  direction  and  are  on  the  X  -X.- 

plane.   Each  boom  has  length  L  and  rectangle-tubular  cross 
section  of  flexural  rigidities  EI.  about  axis  X.  and  EI?  about 

axis  X   (for  C,D)  or  about  X   (for  E,  F) . 

B  has  mass  M  and  the  mass  distribution  in  B  is  assumed  to  be 
axisymmetric  about  X. .   Therefore  among  the  principal  moments  of 

inertia  B  ,  B  ,  B   of  B  for  the  mass  center  B*  associated  with 

X  ,  X  ,  X_ ,  respectively,  B   and  B_  are  the  same.   Each  boom  has 


mass  m,  and  the  mass  of  a  magnetometer  attached  at  the  tip  of 
each  boom  is  m. 

The  orientation  of  B,  or  the  X.,  X2,  X3,  in  a  Newtonian 

reference  frame  N  can  be  specified  in  terms  of  three  external 
coordinates,  e    e  ,  8  .   Let  e^  e2,  &3    be  Euler  angles  such 

that  X.  ,  X0 ,  X^  coordinates  system  is  brought  into  an  orientation 

in  N  by  successive  rotations  a,  ^,  e  3  about  X.,  X  ,  X  from  the 

N  ,  N2,  N3  direction. 

The  deformed  configuration  of  the  satellite  is  described 
using  4n  internal  coordinates  q  ,  q  ,  ...  q4n,  in  addition  to  the 

external  coordinates.   The  general  deformation  of  a  boom  C  is 
described  using  the  n  internal  coordinates. 


n 

Wc(x,t)  =  ^W^t)  (1) 


i=l 

where    W     (x,t)     is    the    deflection    of    a    boom    cross-section    at 

position   x    from   the   boom   root    in    X      X      X^    coordinates,     and  $± 

i  =  l,2,...n  are  admissible  functions  chosen  to  be  the  same  as 
the  normal  modes  of  the  boom  as  a  cantilever  beam  with  an  end 
mass  in  this  analysis.  The  general  deformation  of  booms  D,  E,  F, 
are  described  using  the   same    form   of    equation    (1)    with    q.  ,     i    = 

n+l,...,2n    for  boom   D,    with   q.,    i  =   2n+l,...,3n   for  boom  E,    with 

q. ,    i   =   3n+l,...,4n   for  boom  F,    respectively. 

It   is   assumed  that  no   external    forces   and  moments   are 
applied.      Due  to  the   internal    forces   of    interaction   between    the 
particles    of    systems,    it   is  presumed  that   energy  dissipation  can 
be    described    by    functions    D     (q,q)     which    are    functions    of 

qx / • • • / q4n /    ^i ' • • • ' §4n     sucn  tnat 

Dr(q,0)    =0,  r  =    1,2, ...,4n  (2) 

Dr(q,q)     qr      ^    0  (3) 

where    the    equality    sign    of    equation     (3)     holds     only    when 
q.,...,q        all   vanish. 


1-2   Stability  Conditions 

In  torque-free  motion,  the  satellite  system  S  can  move  such 
that  the  angular  velocity  of  B  in  the  Newtonian  reference  frame  N 
has  a  constant  magnitude  q      and  is  parallel  to  X.  which  remains 

fixed  in  N  while  the  booms  remain  normal  to  X..   This  motion  is 

defined  as  a  simple  spin. 

A  simple  spin  is  called  stable  if  the  system  resumes  this 
motion,  or  acquires  a  motion  closely  resembling  it,  subsequent  to 
every  sufficiently  small  disturbance  that  leaves  the  angular 
momentum  unaltered.  If  the  system  returns  to  the  original  simple 
spin  state  after  a  disturbance  as  time  t  approaches  infinity, 
then  the  simple  spin  is  said  to  be  asymptotically  stable. 

Stability  conditions  can  be  established  by  stating 
conditions  under  which  the  function  F(e,q,q)  defined  as 

F  -   2(E-EQ)  (4) 

is  a  Liapunov  function  (Ref .  7)  ,  where  E  is  the  total  energy  of 
the  system  and  E   is  the  total  energy  of  the  system  during  the 

simple  spin.   It  was  shown  in  Ref.  4  that  for  a  simple  spin  of  an 
axisymmetric  system, 


2    T°   T° 

F  =  —  [Ji±  (  Q22+  032  )  }  +  Z°,re  qrqs  +  M°  qrqs  +  O3(0,  q,  q)  (5) 

1°        1° 

s      t 


where  H  is  the  magnitude  of  angular  momentum,  I   and  I.  are  the 

principal  moment  of  inertia  of  the  system  about  the  axial  axis 
and  about  the  transverse  axis,  respectively,  0_(e,q,q)  indicates 

terms  of  third  or  higher  degree  in  its  arguments,  and  superscript 
o  means  its  value  at  the  simple  spin  state. 

M    is  an  element  of  an  inertia  like  matrix  defined  as 
rs 

M  =  <mp,  .  p,  >  -  <mp  x  p,  >  .  f  -  <mp  x  p,s>  (6) 

r,  s  =  1,2, ,4n 


8 


where  the  symbol  <  >  denote  summation  over  all  particles  of  the 
system  S,  p  is  the  position  vector  of  a  particle  relative  to  the 
mass  center  of  F,  p,   indicates  its  derivative  about  q  ,  m  is  the 

—  1   .  . 

mass  of  the  particle,  and  I     is  the  inverse  of  the  inertia 

dyadic  I. 

Z.    is  a  second  derivative  of  an  energylike  function  Z(q) 
rs 

with  respect  to  q  and  q  ,  which  is  defined  as 

Z(q)  =  V(q)+iH2r1(q)  (7) 

where  V(q)    is   the  potential   energy   function   of   S. 

From    the    Liapunov    stability    theory,     the     stability 
conditions   are 

F      <         0  (8) 

with  equality  sign  holds  only  when  q  ,  ...,q   all  vanish, 


1°   -L°    >   0  (9) 


s      ~t 


and 


Z,     is  positive  definite.  (10) 

27S 


The  first  condition,  equation  (8)  is  satisfied  from  the 

energy  dissipation,  equations  (2)  and  (3) .   The  second  condition, 

equation  (9),  is  the  well  known  "Maximum  Axis  Rule"  of  the 

stability  condition  of  a  rigid  body,  and  can  be  called  the 

"external  stability  condition"  of  a  deformable  body  system.   This 

condition  can  be  satisfied  when  the  boom  length  L  is  longer  than 

the  critical  length  L    ,  which  is  about  60  inches  for  the  NPS 

mm 

mini-satellite  (ORION) . 

The  third  stability  condition,  equation  (10) ,  is  due  to  the 


flexible  deformation  of  booms  and  can  be  called  the  "internal 
stability  condition".  It  is  shown  in  Ref.  4  that  for  a  stress- 
free  spin  about  axis  X.,  that  is, 

V°,r  =  0      ;    r  =  1,2, ,  4n  (11) 

or 

I°rr  =  0      ;    r=l,  2, ,  4n  (12) 


from  the  associated  eqations  governing  the  internal  coordinates, 
Z°     can  be  expressed  as  follows  with  H   replaced  by  I.,  fl 

Z,rs  =  V°,rs  -  Q2{  J-  (fnj  f12s  +  fu  r  fnJ  +  jtn„)  ;     r  =  l,2,...4n  (13) 

S    t 


The  simple  spin  motion  of  the  NPS  satelite  system  S  is  a 
stress-free  spin,  since  it  satisfies  the  condition,  equation  (11) 
or  equation  (12),  by  choosing  the  internal  coordinates  as  in 
equation  (1)  .   Moreover,  by  this  choice  of  internal  coordinates 
using  the  normal  vibration  modes  of  the  cantilevered  boom,  the 

components  of  Z ,    in  equation  (13)  can  be  easily  calculated, 

using  the  orthonormal  relations  between  normal  modes  of  the  boom: 


m<t>r(L)<t>s(L)  +  J<j>r(x)(|>s(x)  n  dx  =  5^  (14) 


Jekv 


(x)(J)s(x)  dx  =  5  k  (15) 


10 


where  6    is  the  Kronecker  delta  of  r,  s  =  l,2,...,n,  and  m 
and  k    r  =  l,2,...,n,  are  modal  masses  and  modal  stiffnesses  of 

J.  J-  a 

the  vibration  modes,  y  is  the  mass  per  unit  boom  length. 


V°.    =k  5    +IL, 

rs         rr  rs  G 


R+L 

=  nVD?8„  +  Q2)j  (m+(R+L-x)}  ^(X)  <^(x)  dx 


(16) 


where    K_    is    the   geometric    stiffness    of    the    boom    due    to    the 

centrifugal    force   of  the   spin,    and   w      is   the   natural    frequency  of 

the   r  mode. 

R+L 

*I2.r  =  JX  VX)  V  dx  +  L  ♦rCL)  m  ( 17  ) 

R 

r  =  vibration  modes   in  X.-   X?   plane   for  booms   C   and   D 

R+L 

^3,8  =  J  X  ^  V   dX  +  L  V^  m  (  18  ) 

R 
s  =  vibration  modes  in  X.-  X  plane  for  booms  E  and  F 

I1t  =  m  8  fTa\ 

11, rs    rr  rs  (19) 


r  =  vibration  modes  in  X  -  X3  plane  for  booms  C,  D,  E, 
and  F 


,o 


When  all  ingredients  required  for  the  evaluation  of  Z,    are 

available  using  equations  (13)  ,  (16-19)  ,  the  internal  stability- 
conditions  can  be  examinded  as  followes: 


11 


Define       A.  ,     i  =  1,2, 4n>      as 


M?n 


(20/a) 


A        A 

A2  = 


A       A 

A3  = 


z° 

^  '11 


z°, 


21 


^  '11 


z 


21 


z° 

^  'i: 


Z° 


22 


z° 


z°, 


22 


z°  z° 


32 


z°, 
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z° 


23 


z°, 


33 


and  so  on. 

The   internal    stability  condition   is 


(20, b) 


(20/C) 


A.  >0 

i 


for  all     i  =  1,2,3, ,4n 


(21) 


II-3      Stability   of  the  NPS  Mini-Satellite    (ORION) 

The   stability   conditions   of   the   simple   spinning  motion   of 
the  NPS   mini-satellite  between    the    angular   velocity  n       and    the 
length    of     a    boom    L     for    other    parameters    defined    in    the 
preliminary  design    (See  Ref .     1)  .       The    numerical    value    of    these 
parameters   are   shown   in  Table    1. 
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Table  1  -  Parameters  of  NPS  mini-satellite 

Radius  of  B  R  =  9.5  inch 

Mass  of  B  M  =  7.542  slug 

Axial  moment  of  inertia  of  B  B1  =  2 . 18  slug-ft 

Transverse  moment  of  inertia  of  B  B2  =  B3  =  6.78  slug-ft 

Area  of  a  boom  cross-section  A  =  0.212  in 

Area  moment  of  inertia  of  a  boom 

about  X  -  X  plane  I1  =  0.01668  in 

Area  moment  of  inertia  of  a  boom 

about  X.-  X2  plane  I2  =  0.06328  in 

Young's  modulus  of  the  boom  E  =  10.3  X  10   Psi 

Mass  of  a  boom  per  unit  length  v  =  0.000415  (Slug/ in) 

Mass  of  a  magnetometer  at  the  end 

of  a  boom  m  =  0.0621  slug 


For  the  design  value  of  L  =  78.5  inches  the  cantilever  beam 
model  of  a  boom  with  an  end  mass  was  analyzed  using  the 
MSC/NASTRAN  code  for  the  natural  frequencies  and  mode  shapes. 
The  natural  frequencies  are  shown  in  Table  2,  and  the  finite 
element  model  and  the  first  3  mode  shapes  are  shown  in  Figures  4 
to  7. 
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Table  2  -  Natural  Frequencies   of  a  Boom 


Mode 


Natural  Frequency  (Hz) 


Mode  Shape 


1 
2 
3 
4 
5 
6 


2.104 

4.097 

23.426 

45.569 

73.656 

142.980 


1-2  plane  bending 
2-3  plane  bending 
1-2  plane  bending 
2-3  plane  bending 
1-2  plane  bending 
2-3  plane  bending 


When  the  angular  velocity  of  simple  spin  ft  is  increased 
from  very  small  value,  it  can  be  seen  from  stability  conditions, 
equations  (20-21)  and  from  equations  (13),  (16-19)  that  the 
stability  condition  will  be  violated  first  when  ft  approaches  the 
first  natural  frequency  &         since  for  higher  modes  u  ?,  i  = 

2,3,...,  are  much  greater  than  ri?  and  the  matrix  Z,    become 

diagonally  dominant.   Therefore,  the  stability  conditions  can  be 
examined  simply  by  reducing  the  number  of  internal  coordinates  q. 

to  the  first  mode  of  each  boom,  q  ,  q  ,  q  ,  q  . 

Inserting  the  expressions  of  its  components,  equations  (16) 
-  (19)  into  equation  (13) ,  the  matrix  Z,    becomes  the  following 

£  fi) 

form: 


Z°,    = 

rs 


co2-  aQ2 

i 

0 

-(to2 

0 

co2  -  aQ.2 

l 

0 

-po2 

0 

co2  -  a£22 

l 

0 


-[to2 


0 

■pn2 

0 


co2  -  a£22 

i 


(22) 
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where 


R+L  r+l 


a  =  — ^  {J  x  (^(x)  p.  dx  +  L  (^(L)  m}2-  Ji  (m+u.(R+L-x)}  {(^(x^dx         (23) 


r-i:  R  R 


s     t 


and  R+L 


P  =  —  {  [ x  4).  (x)  |i  dx  +  L  0  (L)  m  } : 

To  To     J         l  1 


(24) 


r-L  R 


s      t 


with    <L  (x)    is   the  mode   shape   of  the   first  mode. 


The   stability  conditions,    equations    (20)    -    (21) ,    become 

(0?  -ctf22>0  (25,a) 


(of-afi¥>0  (25/b) 


(co2  -  aQ2){(co2  -  aft2)2  -  (32  Q4}  >  0  (25, c) 


{(co2-aQ2)2-(52Q4}2>0  (^  ^ 

1  (25,0) 


Therefore,    the    internal    stability    conditions    is   reduced   to   the 
condition 


co2 
<  - 


tf-  -1 


tt+P  (26) 


from  the   condition    (25, c). 
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As  the  length  of  a  boom  L  increases,  the  first  natural 
frequency  w  decreases  as  follows; 


3EI 


(m+0.23uL)L3 


(27) 


The  denominator  in  equations  (23) ,  (24) ,1-1.    is 

2     2 

proportional  to  (L   -  L  .  )  and  increases  from  zero  for  L  >  L  .  , 

c  mm  mm 

which  is  the  external  stability  condition.  Therefore,  the 
stability  boundary  in  the  stability  diagram  varies  as  shown  in 
Figure  8.  For  the  designed  value  of  L  =  78.5  inches  the 
stability   condition   is 

Q<  8.658  (rad/sec)  (28) 

However,  it  was  assumed  that  each  boom  is  a  single  piece  and 
rigidly  connected  to  the  rigid  main  body  B.  Since  each  boom 
consist  of  three  beam  and  the  connection  between  the  boom  and  the 
body  B  is  not  ideally  rigid  the  critical  angular  velocity  ft 

may  be  much  lower  than  the  value  calculated,  equation  (28) , 
depending  on  the  design  of  joints  of  booms. 
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III.   Conclusions 

The  stability  conditions  of  a  simple  spin  motion  of  the  NPS 
mini-satellite  (ORION)  are  derived  analytically  using  Liapunov 
stability  theory  and  by  choosing  the  natural  vibration  modes  as 
admissible  functions  of  internal  coordinates  of  the  system.  By 
applying  this  conditions  to  the  preliminary  design  of  the 
satellite  it  is  shown  that  the  angular  velocity  .of  the  simple 
spin  should  be  limited  to  the  value  of  8.685  rad/sec  for  the  boom 
length  L  =  78.5  inches.  If  the  boom  length  is  changed  the 
critical  angular  velocity  Q      also  changes  and  the  stability 

conditions  should  be  re-examined  using  the  results  derived  in 
this  study. 

It  should  be  recognized  that  the  stability  conditions  were 
derived  for  the  simplified  ideal  model  under  the  following 
assumptions;  1)  the  mass  distribution  of  the  cylindrical  main 
body  is  axisymmetric;  2)  the  connections  between  the  booms  and 
the  main  body  and  the  main  body  itself  are  rigid;  3)  each  boom 
acts  as  a  single  piece  and  the  locking  mechanism  of  the  boom 
should  be  ideal. 

Due  to  the  highly  nonlinear  behavior  of  a  three  dimensional 
motion  of  a  flexible  body,  the  stability  analysis  is  valid  for 
only  small  disturbances.  Therefore,  stability  of  dynamic  motion 
of  a  satellite  for  general  environments  should  be  investigated  by 
dynamic  simulation  using  computers  and  proper  dynamic  models. 
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Figure  1.   The  NPS  Mini-Satellite  (ORION) 
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Figure  2.   Random  Vibration  Level  in  GAS  Canisters 
(NASA   G.  A.  S.  Manual,  1984,  P.  57) 
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Figure  3.   Dynamic  Model  of  The  NPS  Satellite 

During  a  Simple  Spin  Motion  in  Space 
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Figure  4 .   The  Finite  Element  Model  of  a  Boom 

(X]_,  X2  ,  X3   show  the  Orientation  of  the  Satellite) 
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Figure  5.   The  Mode  Shape  of  the  First  Mode  (f  =  2.104  Hz) 
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Figure  6.   The  Second  Mode  (f  =  4.097  Hz) 
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Figure  7.   The  Third  Mode  (f  =  2  3.43  Hz) 
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Figure  8.   Stability  Diagram  of  The  NPS  Mini-Satellite 
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